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Several  anisotropic  models  for  the  three-dimensional  spectra  of  velocity 
fluctuations  in  shear-dominated,  atmospheric  surface-layer  turbulence  are 
presented  and  discussed.  These  models  include  a  generalized  Gaussian 
model,  two  forms  of  generalized  von  Karman  models;  the  Kristensen, 
Lenschow,  Kirkegaard,  and  Courtney  model;  and  the  Mann  model.  The  dis¬ 
tinction  between  a  "top-down"  and  "bottom-up"  approach  to  the  design  of 
a  model  is  discussed,  and  how  the  bottom-up  approach  generally  leads  to 
more  satisfactory  models  is  shown.  The  effects  of  turbulent  anisotropy  on 
acoustic  propagation  are  explored  by  calculating  mutual  coherence  func¬ 
tions  (MCFs)  (describing  the  coherence  of  a  propagating  acoustic  wave)  for 
the  different  models.  Anisotropy  effects  have  been  found  that  they  can  be 
quite  significant,  even  when  the  separation  between  the  acoustic  sensors  is 
small. 
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1.  Introduction 


Acoustical  arrays  can  be  used  to  determine  the  direction  of  a  source.  The 
principle  is  simple:  through  the  use  of  phase  differences  across  the  array, 
the  orientation  of  the  incident  wave  fronts  can  be  determined.  If  little  re¬ 
fraction  exists,  the  normals  to  the  wave  fronts  correspond  to  the  direction 
of  the  source.  An  important  limiting  factor  in  this  process  is  the  distortion 
of  the  propagating  wave  fronts  caused  by  atmospheric  turbulence.  The  tur¬ 
bulence  introduces  random  perturbations  to  the  wave  fronts,  so  that  they 
are  no  longer  perfectly  spherical  waves  radiating  outward  from  the  source 
(fig.  1).  In  such  situations  that  the  signal-to-noise  ratio  (SNR)  is  high,  the 
turbulence  effect  is  often  the  dominant  factor  affecting  performance  of  the 
array  (Wilson,  1997a). 

The  turbulence  effect  can  be  quantified  by  calculating  the  mutual  coher¬ 
ence  function  (MCF),  which  describes  the  coherence  between  signals  at  a 
pair  of  sensors  as  a  function  of  their  spatial  separation.  The  MCF  for  prop¬ 
agation  through  small-scale  (inertial  subrange)  turbulence  was  determined 
long  ago  and  has  been  applied  successfully  to  optical  and  other  types  of 
wave  propagation  (see,  for  example,  Tatarskii  (1971)).  Treatment  of  the  in¬ 
ertial  subrange  is  reasonably  straight  forward,  primarily  because  the  statis¬ 
tics  of  the  small  eddies  are  locally  isotropic  (independent  of  orientation)  and 
homogeneous  (independent  of  position).  In  the  case  of  low-frequency  acous¬ 
tical  propagation,  however,  large-scale  (energy  subrange)  turbulence  plays 
an  important  role.  Hence  it  is  necessary  to  account  for  such  intrinsic  fea¬ 
tures  of  the  energy  subrange  as  anisotropy  and  inhomogeneity.  An  example 


Figure  1.  Random  distortions  in  acoustic  wave  fronts  caused  by  atmospheric  turbulence. 
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of  anisotropy  is  that  the  variance  for  turbulent  velocity  fluctuations  in  the 
direction  of  the  mean  wind  is  several  times  larger  than  the  variance  for  ver¬ 
tical  velocity  fluctuations.  An  example  of  inhomogeneity  is  the  dependence 
of  the  variances  on  the  height  above  the  ground.  Similar  manifestations  of 
anisotropy  and  inhomogeneity  are  evident  in  the  turbulent  length  scales. 

In  a  previous  report  (Wilson,  1997b),  I  discussed  a  model  for  inhomo¬ 
geneities  in  atmospheric  turbulence  that  are  introduced  by  "blocking"  at 
the  ground.  This  report  is  intended  as  a  complement  to  the  previous  one.  Its 
focus  is  on  anisotropic  models  for  atmospheric  turbulence  and  on  the  impli¬ 
cations  for  acoustic  wave  coherence  as  quantified  by  the  MCE  As  in  the  pre¬ 
vious  report,  the  turbulence  models  discussed  are  fully  three-dimensional 
(3D),  second-order  statistical  models.  By  this  I  mean  that  the  models  can 
be  applied  to  any  second-order  statistic  of  the  turbulence,  with  the  orien¬ 
tations  of  the  velocity  components,  as  well  as  the  direction  of  the  spatial 
displacement  between  the  two  observation  points,  being  arbitrary. 

Although  the  turbulence  models  discussed  in  this  report  can  be  applied,  in 
principle,  to  a  broad  range  of  atmospheric  conditions,  the  selection  of  the 
model  parameters  described  is  for  the  specific  case  of  a  shear-driven,  neu¬ 
tral,  atmospheric  surface  layer.  Neutral,  in  the  terminology  of  geophysical 
fluid  dynamics,  means  an  absence  of  buoyant  instabilities.  Hence  radiative 
heating  and  cooling  of  the  ground  (and  the  resulting  temperature  gradients 
in  the  overlying  air)  are  neglected.  This  leaves  wind  shear  as  the  dominant 
mechanism  for  generating  the  turbulence. 

One  of  the  points  that  is  demonstrated  in  this  report  is  that  the  modeling 
procedure  becomes  very  important  for  anisotropic  turbulence.  If  the  proce¬ 
dure  is  ill-conceived,  the  result  will  often  be  a  poor  model  that  does  not 
agree  well  with  data  and/or  has  parameters  that  are  nearly  impossible  to 
determine.  There  are  basically  two  possible  procedures,  which  I  call  the 
"top-down  approach"  and  the  "bottom-up  approach."*  These  approaches 
are  illustrated  in  figure  2.  The  top-down  approach  has  been  used  in  the  past 
to  derive  the  popular  isotropic  von  Karman  model.  It  begins  by  postulat¬ 
ing  an  equation  for  the  specific  turbulent  kinetic  energy  (TKE)  spectrum"^  (just 
energy  spectrum  for  short)  of  the  turbulence,  as  was  done  in  the  previous 
report  (Wilson,  1997b): 


^  4r(17/6) 

^  ’  V^r(l/3)(l  +  P£2)17/6- 


(1) 


‘This  usage  of  "top-down"  and  "bottom-up"  should  not  be  confused  with  the  terminol¬ 
ogy  coined  for  atmospheric  boundary  layer  diffusion  by  Wyngaard  and  Brost  (1984). 
^Specific  means  per  unit  mass. 
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Figure  2.  Two  possible  procedures  for  deriving  statistical  turbulence  models. 


The  energy  spectral  equation  has  two  unknown  parameters,  a  length  scale 
i  and  a  variance  (Other  symbols  and  notation  used  in  equation  (1)  are 
discussed  in  sect.  2  of  this  report.)  When  homogeneity  and  isotropy  are 
assumed,  a  well-known  set  of  relationships  allows  derivation  of  all  other 
second-order  statistical  quantities,  such  as  correlations  and  spectra,  from 
the  energy  spectrum.  For  example,  the  one-dimensional  (ID)  spectrum, 
Fl  (k),  which  is  measured  by  a  stationary  sensor  when  the  fluid  advects 
the  turbulence  past  the  sensor,  is  related  to  the  energy  spectrum  through 
the  equation  (Batchelor,  1953) 


E{k) 


ri  dFL  {ky 

dk  k  dk 


It  can  be  easily  verified  that 


^  4r(5/6)  aH 

>  V^r(l/3)(i  +  p£2)5/6' 


(2) 


(3) 


The  key  to  determining  the  unknown  parameters  {£  and  in  this  case)  is  to 
derive  one  or  more  equations  for  observable  quantities,  such  as  Fl  (k).  The 
parameters  can  then  be  deduced  by  fitting  the  derived  equations  to  obser¬ 
vations.  The  reason  for  calling  this  approach  top-down  is  that  the  starting 
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point  for  the  modeling  procedure  is  the  fundamental,  underlying  function: 
the  energy  spectrum. 

The  starting  point  in  the  second,  or  bottom-up,  approach  is  the  observed 
data.  An  equation  for  the  ID  spectrum  is  then  postulated,  and  its  parame¬ 
ters  are  chosen  according  to  the  data.  Then  the  more  fundamental  function, 
the  energy  spectrum,  is  determined  in  some  manner  from  the  ID  spectrum. 
This  is  the  approach  that  has  usually  been  adopted  in  the  past  to  derive  the 
popular,  isotropic  Gaussian  model. 

In  the  comparison  of  the  top-down  and  bottom-up  approaches,  it  is  impor¬ 
tant  to  understand  that  the  energy  spectrum  cannot  be  determined  from 
the  ID  spectrum  without  making  some  set  of  assumptions,  such  as  homo¬ 
geneity  and  isotropy,  and  then  deriving  a  relationship  such  as  equation  (2). 
The  process  of  determining  the  ID  spectrum  from  the  energy  spectrum  is 
generally  much  more  straight  forward,  involving  integrations  that  are  valid 
regardless  of  the  turbulence  properties. 

For  homogeneous  and  isotropic  turbulence,  the  relationship  between  the 
energy  spectrum  and  the  ID  spectrum  is  so  simple  (eq  (2)),  it  does  not  re¬ 
ally  matter  whether  one  uses  a  top-down  or  bottom-up  approach;  it  is  easy 
to  convert  back  and  forth  between  the  energy  spectrum  and  the  ID  spec¬ 
trum.  When  anisotropy  is  considered,  however,  this  conversion  becomes 
much  more  difficult.  Therefore  some  thought  must  be  given  to  the  distinc¬ 
tion  between  the  top-down  and  bottom-up  approaches  in  the  anisotropic 
case.  We  see,  in  fact,  that  although  the  top-down  approach  is  more  conve¬ 
nient  from  a  mathematical  standpoint,  it  is  very  difficult  to  derive  a  model 
agreeing  well  with  observations  when  this  approach  is  used.  On  the  other 
hand,  the  bottom-up  approach  becomes  mathematically  cumbersome,  al¬ 
though  Kristensen  et  al  (1989)  have  devised  a  systematic  process  for  the 
derivation.  The  main  reason  for  adopting  the  bottom-up  approach  is  that  it 
allows  for  creation  of  a  model  having  good  agreement  with  observations. 

I  consider  six  anisotropic  models  in  this  report.  The  reader  may  wonder 
why  it  is  necessary  to  consider  such  a  multitude.  I  am  not  suggesting  that 
all  of  the  models  are,  in  the  end,  useful  ones.  Rather  I  introduce  the  large 
number  of  models  in  this  report  to  document  how  anisotropic  models  can 
be  developed  from  various  starting  points,  and  to  compare  their  relative 
merits  and  fidelity  to  atmospheric  measurements.  By  the  time  we  finish, 
it  should  be  clearer  which  models  are  actually  worth  retaining  for  future 
studies. 

The  first  model  I  consider,  a  set  of  empirical  equations  owing  to  Kaimal  et 
al  (1972)  (sect.  3),  is  actually  valid  only  for  ID  spectra,  with  the  direction  of 
the  displacement  along  the  axis  of  the  mean  wind.  In  this  report,  our  main 
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use  for  Kaimal  et  al's  spectral  equations  is  as  a  reference  point  for  choosing 
parameters  needed  by  the  other,  complete  models  derived  in  sections  4  to 
7.  For  better  or  worse,  the  Kaimal  et  al  equations  are  used  as  the  "ground 
truth." 

The  second  model  I  consider  (sect.  4),  which  is  original  with  this  report, 
is  essentially  the  top-down  generalization  of  the  popular  anisotropic  von 
Karman  model.  The  third  model  (sect.  5),  also  original  with  this  report, 
is  a  bottom-up  generalization  of  the  anisotropic  von  Karman  model.  The 
fourth  model  I  consider  (sect.  6)  is  due  to  Kristensen  et  al  (1989).  The  Kris- 
tensen  et  al  model  is  actually  quite  similar  to  the  bottom-up  anisotropic  von 
Karman  one,  although  it  is  somewhat  more  complicated  and  flexible.  The 
fifth  model  I  consider  is  an  anisotropic,  bottom-up,  generalization  of  the 
popular  Gaussian  model  (sect.  7).  The  last  model,  the  Mann  (1994)  model 
(sect.  8),  is  entirely  different  from  the  previous  ones:  it  is  a  first-principle 
treatment  for  sheared  turbulence.  The  Mann  model  is  rather  complex,  but  it 
contains  certain  realistic  features  of  sheared  turbulence  that  the  other  mod¬ 
els  do  not  capture.  In  section  9  of  this  report,  example  calculations  of  MCFs 
and  source  direction-finding  performance  for  acoustical  arrays  are  made 
using  the  various  anisotropic  turbulence  models. 
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2.  Preliminaries 


This  section  reviews  many  of  the  basic  correlation  and  spectral  relation¬ 
ships  involved  in  modeling  of  turbulence  statistics  and  defines  the  nota¬ 
tions  used  in  subsequent  sections  of  the  report.  It  is  assumed  that  the  reader 
has  a  basic  understanding  of  statistics  and  Fourier  transformation. 

The  correlation  function  tensor  is  defined  as 

Rij  {x,x  +  r)  =  {ui{x.)uj  (x  +  r)) ,  (4) 

where  the  angle  brackets  indicate  ensemble  averaging,  and  Ui  and  Uj  are 
velocity  fluctuations  in  the  fth  and  jth  directions.  In  this  report,  I  follow 
the  usual  convention  for  orientation  of  the  coordinates  in  micrometeorol¬ 
ogy:  the  xi-axis  is  aligned  with  the  direction  of  the  mean  wind  (assumed 
horizontal),  the  X2-axis  is  the  horizontal  axis  perpendicular  to  the  mean 
wind,  and  the  xa-axis  is  vertical. 

In  homogeneous  turbulence  (an  assumption  made  throughout  this  report), 
the  correlations  depend  only  on  the  separation  between  the  observation 
points: 

Rij  (r)  =  Rij  (x,  X  -h  r) .  (5) 

The  correlation  function  tensor  forms  a  Fourier-transform  pair  with  the  3D 
spectral  density  tensor  (k): 

i^)  =  J  (k)  exp  (ik  •  r)  d^k,  (6) 

(k)  =  7-^  /  Rij  (r)  exp  (-ik  •  r)  d^r.  (7) 

yZlT)  J 

(The  notation  f  d^k  in  eq  (6)  means  that  the  integration  is  to  be  carried  out 
over  the  entire  3D  k-space,  and  likewise  r-space  for  eq  (7).) 

The  covariance  tensor  is  defined  as 

afj  =  (ui  (x)  Uj  (x))  =  Rij  (0)  .  (8) 

By  evaluating  equation  (6)  for  r  =  0,  we  see  that  the  covariance  can  also  be 
found  from  the  3D  spectral  density  by  integrating  over  the  wave  number 
domain: 

4.  =  J  CD,,  (k)  d^k.  (9) 
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Normally  in  experiments  the  ID  spectral  density  tensor  Qij  (k)  is  measured. 
This  quantity  is  defined  as 

1 

Qij,n  {k)  =  —  Rij  (re„)  exp  {-ikr)  dr,  (10) 

J  —  oo 


where  e„  is  the  unit  vector  along  the  Xn-axis.  In  the  case  of  atmospheric 
turbulence,  the  measurement  usually  is  performed  from  a  stationary  tower, 
and  Taylor's  hypothesis  is  invoked  to  convert  the  time  series  to  a  spatial 
one  (Panofsky  and  Dutton,  1984).  In  this  case,  e„  points  in  the  direction  of 
the  mean  wind. 

Alternatively,  we  can  determine  the  ID  spectral  density  from  the  3D  spec¬ 
tral  density.  Suppose  we  apply  a  ID  Fourier  transform  to  equation  (6)  with 
respect  to  the  Xn-axis  and  then  evaluate  the  result  with  r  =  re„.  We  find 


1 


Rij  (re„)exp(-f/cnr) 


dr 


^ij  (k)  d^k±, 


(11) 


where  k  =  kriGn  +  k^,  k^  =  i  -|-  k±^2^^^2,  and  the  subscripts  "T,  m" 

indicate  the  directions  perpendicular  to  the  Xn-axis.  The  left  side  of  this 
equation  is  the  same  as  the  right  side  of  equation  (10),  and  so  we  have  the 
result 

0p>  {k)  =  J  {ken  -h  k_L)  d'^k±.  (12) 

For  brevity,  I  often  use  Fi  (k),  Ft  (k),  and  Fy  (k)  in  this  report,  instead  of 
©11,1  (k),  022,1  (k),  and  ©33, i  (k).  This  set  of  ID  autospectra,  having  the 
wave  number  vector  aligned  with  the  direction  of  the  mean  wind,  is  the 
one  observable  using  a  stationary  sensor. 

An  important  quantity  in  wave  propagation  is  the  two-dimensional  (2D)  cor¬ 
relation  function  bij^n  (r±): 

I  roo 

kij,n  ('^x)  —  X  /  ki'i'j  {^n^n  T  '^x)  dCn,  (13) 

J —  oo 


where  =  rx,iex,i  -l-  xx,2ex,2-  In  wave  propagation  problems,  the  Xn- 
axis  is  the  direction  of  propagation,  and  usually  only  the  wind  components 
parallel  to  the  direction  of  propagation  are  relevant.  Hence  i  =  j  =  n  is  the 
main  interest.  A  relationship  between  the  2D  correlation  and  3D  spectral 
density  can  also  be  derived,  by  taking  the  ID  Fourier  transform  of  equa¬ 
tion  (6)  with  respect  to  the  x^-axis,  and  then  evaluating  the  result  with 
kn  =  0.  We  find 


(rx)=/' 


(k_L)  exp  (ik_L  •  r_L)  d‘^k±. 


(14) 
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In  the  context  of  wave  propagation,  it  is  also  useful  to  define  a  2D  structure 
function  dij^n  (r±)  as 

dij,n  (r_L)  =  2  [hij^n  (0)  -  (r±)]  •  (15) 


The  2D  structure  function  is  the  2D  counterpart  of  the  usual  (3D)  structure 
function  encountered  in  turbulence  theory. 

We  define  the  integral  length  scale  as 

I  roo 

d-‘ij,n  —  2~  /  (^®«)  df.  (lb) 

Jo 


Various  useful  relationships  between  the  integral  length  scale,  the  spectral 
densities,  and  the  2D  correlation  function  can  be  derived  from  the  Fourier 
transform  relations.  First  note  that  when  i  =  j,  the  correlation  function  is 
an  even  function  and  equation  (16)  can  be  rewritten 


d-‘ii,n  — 


1 

— ^  /  Rii  (re„)  dr. 

J-oo 


By  taking  the  ID  Fourier  transform  of  equation  (6)  with  respect  to  the  Xn- 
axis,  and  then  evaluating  the  result  with  kn  =  0  and  =  0,  it  can  be  shown 


that 


r-  -  — 

^ii.n  —  9 


<^ii  (k_L)  d‘^k±  . 


(17) 


By  setting  /c  =  0  in  equation  (12),  we  therefore  have 


L, 


=  4©^ 

erf, 


(0) 


(18) 


Furthermore,  by  setting  =  0  in  equation  (14),  we  have 


TT 

Lii,n  —  2  (0) 


(19) 
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3.  Kaimal  et  al  Model 


Kaimal  et  al  (1972)  developed  a  simple  set  of  empirical  formulas  for  the 
ID  velocity  spectra.  Their  equations  are  actually  for  the  frequency,  rather 
than  the  wave  number,  spectra,  such  as  discussed  in  section  2.  With  the  as¬ 
sumption  of  Taylor's  hypothesis  though  (which  states  that  the  turbulence 
is  "frozen"  as  it  moves  past  the  sensor  at  the  mean  wind  speed  U),  their  fre¬ 
quency  n  is  simply  equal  to  Uk/2TT.  As  a  result  Kaimal  et  al's  equations  21a 
to  21b  for  the  spectra  of  the  horizontal  velocity  components  become,  in  the 
notation  of  this  report. 


Fl  {k)  =  011,1  (fc) 


_ bizu^ _ 

Att  (1  +  ci\k\  ’ 


(20) 


Ft  {k)  =  022,1  {k)  = 


Att  (1  -|-  C2  \k\  z/27r)^'^^ 


where  bi  =  105,  ci  =  33,  62  =  17,  and  C2  =  9.5.  The  equation  for  the  vertical 
velocity  spectrum  is 


Fv  (fc)  =  033,1  {k)  = 


Att  1  -|-  C3  I 


where  63  =  2  and  C3  =  5.3.  Kaimal  et  al's  equations  are  plotted  in  figures  3 
to  5. 


Kaimal  et  al's  equations,  since  they  describe  only  the  ID  spectra,  do  not 
constitute  a  complete  3D  model  for  anisotropic  turbulence.  However,  in 
the  remainder  of  this  report,  they  are  used  frequently  as  a  reference  point. 
They  are  compared  to  the  various  complete  models  that  are  developed  fur¬ 
ther  on,  and  the  parameters  in  the  complete  models  (such  as  variances  and 
length  scales)  are  chosen  to  produce  a  good  agreement  with  Kaimal  et  al's 
equations. 

The  variances  corresponding  to  the  Kaimal  et  al  spectra  can  be  found  by 
integrating  the  spectra  from  k  =  —00  to  00  (eq  (1)).  The  result  for  the  hori¬ 
zontal  velocity  components  is 


2  ^biul 
O',-,-  —  — — 
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10° 


Figure  3.  One-dimensional  turbulence  spectra  for  streamwise  (longitudinal)  velocity 
fluctuations. 


Figure  4.  One-dimensional  turbulence  spectra  for  horizontal,  cross-stream  (lateral)  velocity 
fluctuations. 
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For  the  vertical  velocity. 


— 

<^33  — 


5  sin  (27r/5) 

The  computed  values  for  the  variances  are  shown  in  table  1. 

The  integral  length  scales  can  be  determined  easily  using  equation  (18).  We 
find 

bizul 


TmjI  — 


4u£. 


The  resulting  values  for  the  integral  length  scales  are  also  shown  in  table  1. 


Figure  5.  One-dimensional  turbulence  spectra  for  vertical  velocity  fluctuations. 


Table  1.  Values  for  variances  a%  and  sfreamwise  integral  length  scales  La^i  corresponding 
to  Kaimal  model. 


i 

2 

1 

4.77m* 

5.502 

2 

2.68m; 

1.59z 

3 

1.46m; 

0.3432 
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4.  Von  Karman-Based  Anisotropic  Model,  Top-Down  Approach 


In  this  section,  I  attempt  to  develop  a  model  for  anisotropic  turbulence  us¬ 
ing  a  top-down  approach.  The  model  equations  are  based  upon  those  used 
for  the  popular  isotropic  von  Karman  model. 


4.1  Anisotropic  Spectral  Density  Tensor 

Kristensen  et  al  (1989)  have  proposed  the  following  equation  for  the  3D 
spectra  in  anisotropic  turbulence: 

3..,  (k)  =  ^  .4™  (k)  ({„.  -  (i„,  -  ,  (23) 

m=l  ^  ^  ^  ^ 


Their  equation  satisfies  the  incompressibility  condition  for  the  turbulence. 
Note  that  when  Ai  (k)  =  A2  (k)  =  A3  (k)  =  A  (k),  the  equation  reduces  to 
the  familiar  equation  for  isotropic  turbulence  (e.g.,  Batchelor,  1953) 


where 

=  (25) 

and  E  {k)  is  the  energy  spectrum  for  turbulent  kinetic  energy  (TKE). 

An  important  point  regarding  equation  (23)  is  that  it  is  only  an  idealized 
form  for  the  3D  spectra  in  anisotropic,  homogeneous  turbulence.  In  devis¬ 
ing  it,  Kristensen  et  al  neglected  many  of  the  assymetries  that  are  possi¬ 
ble  in  anisotropic  turbulence.  For  example,  equation  (23)  has  the  symmetry 
property  T>i3  (/ci,  k2,  k^)  =  — T>i3  {ki,  k2,  —k^).  Hence  when  T>i3  is  integrated 
from  kz  =  —00  to  00,  the  result  is  zero.  This  implies  that  the  ID  spectrum 
©13,1  {k)  is  zero  for  all  k  (see  eq  (12))  and  therefore  also  that  the  covari¬ 
ance  {uiUs)  is  zero.  But  negative  {uiu^)  is  a  fundamental  characteristic  of  a 
shear  layer,  since  there  is  a  transfer  of  streamwise  momentum  toward  the 
surface.  Despite  such  inherent  shortcomings,  equation  (23)  is  used  to  de¬ 
velop  most  of  the  anisotropic  turbulence  models  in  the  remainder  of  this 
report.  It  is  just  a  first  step  toward  creating  turbulence  models  with  realistic 
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anisotropy;  we  will  see  that  even  this  simple  first  step  causes  very  signifi¬ 
cant  additional  complications.  In  section  8, 1  consider  an  altogether  differ¬ 
ent  model,  the  Mann  model,  that  does  break  away  from  the  constraints  of 
equation  (23),  predicting  a  negative  value  for  (uiUs). 

The  basic  idea  of  the  top-down  approach  is  that  we  build  our  turbu¬ 
lence  model  from  equations  for  the  energy  spectra.  When  the  turbulence 
is  isotropic,  a  flexible  and  useful  form  for  the  energy  spectrum  is 


where 


E{k)  = 


B  [{ij,  +  1)  /2,  u]  (1  +  p£2)(A‘+i)/2-ri^  ’ 


B{x,y)  = 


r (x) r  jy) 
T{x  +  y) 


(26) 


(27) 


is  called  the  beta  function,  and  T  (•)  is  the  gamma  function.  The  reader  can 
verify,  using  equation  (A-1),  that  integration  of  this  energy  spectrum  from 
k  =  0  to  k  =  oo  yields  3a‘^(2  for  all  values  of  p  and  u.  Hence  is  the 
variance  of  a  single  velocity  component.  The  parameter  p  controls  the  slope 
of  the  spectrum  in  the  energy  subrange,  k£  1,  where  E  {k)  oc  /c^.  For  the 
inertial  subrange,  ki^  I,  E  {k)  oc  k~^~‘^''. 

What  is  classically  called  the  von  Kdrmdn  model  corresponds  to  y  =  4,  and 
u  =  1/3.  This  value  for  u  produces  an  inertial-subrange  spectrum  propor¬ 
tional  to  in  agreement  with  Kolmogorov's  scaling  arguments  for  this 

part  of  the  spectrum  (Kolmogorov,  1941).  The  classical  von  Karman  energy 
spectrum  was  given  earlier  as  equation  (1). 

In  this  paper,  I  assume  that  the  Ai  (k)  each  have  y  =  6  and  v  =  1/3.  The 
reason  for  choosing  p  =  6  is  simply  that  this  value  allowed  quantities  such 
as  the  ID  spectral  densities  to  be  derived  in  analytical  form,  whereas  y  =  4 
did  not.  (Actually,  any  value  y  =  2n,  where  n  >  3,  will  work.)  Given  these 
assumptions,  our  equation  for  A*  (k)  is 

(k)  -  3 _ ajkHJ 

4vri?(7/2,l/3)(i  +  P£2)23/6- 

In  this  report,  I  call  this  turbulence  model,  based  entirely  on  equations  (23) 
and  (28),  the  top-down  anisotropic  von  Kdrmdn  model. 

Note  that  six  parameters  are  in  this  version  of  the  anisotropic  von  Karman 
model.  However,  since  the  inertial  subrange  must  be  isotropic,  we  require 
for  k£i  »  1  that 

Ai  {k)  =  A2  {k)  =  As  {k) . 


13 


As  a  result  we  must  have 


(29) 

Let  us  define 

2  ai  -|-  02  + 

^  “  3 

(30) 

and 

,2/3  .2/3  .2/3 

P‘2/3  _  -*^1  ^2  ^  ^3 

3 

(31) 

It  can  then  be  shown  that 

i  =  1,2,3. 

(32) 

This  relationship  will  be  useful  later  in  this  section. 

The  anisotropic  energy  spectra  for  the  top-down  von  Karman  model  are 
plotted  in  figures  6  to  8.  The  parameters  used  for  these  curves  were  chosen 
by  methods  to  be  discussed  in  section  4.5.  As  one  would  expect,  the  A/s 
for  the  top-down  von  Karman  model  have  a  smooth,  single  peak,  reminis¬ 
cent  of  the  isotropic  von  Karman  energy  spectrum.  In  comparison  to  the 
isotropic  von  Karman  model*  (also  shown  in  the  figures),  the  anisotropic 
version  of  Ai  peaks  at  a  higher  wave  number,  whereas  A3  peaks  at  a  lower 
wave  number.  The  peaks  for  A2  roughly  coincide. 
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Figure  6.  Anisotropic  energy  spectra  for  longitudinal  velocity  component. 


‘Equations  for  the  isotropic  von  Karman  model  appear  in  Wilson  (1997a). 
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Figure  7.  Anisotropic  energy  spectra  for  lateral  velocity  component. 


Figure  8.  Anisotropic  energy  spectra  for  vertical  velocity  component. 
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4.2  Variances 


With  the  use  of  equation  (23),  the  spectrum  (k)  is  given  by 
(k)  =  (t)  (^1  -  I  j  (33) 


To  find  the  variance,  we  need  to  calculate  the  integral 

crfi  =  J  4>11  (k)  dk. 


In  spherical  coordinates  {ki  =  k  cos  6  sin  4>,  k2  =  k  sin  6  sin  cj),  k^  =  k  cos  (j)), 
the  integral  becomes 


cjfi  = 


COO  CTV  c27r  r 


/o  Jo  Jo 


A 


(k)  (^1  —  cos^  6  sin^  Jj  +  A2  (k)  (^cos  6  sin  9  sin^ 


+  As  (k)  {cos  6  sin  (p  cos  k^  sin  (p  dO  d(p  dk 


Performing  the  integrations  over  9  and  (p  yields 

/l.'Tr  pOO 

crh  =  —  [8^1  (k)  +  A2  (k)  +  As  (k)]  k"^  dk. 

to  Jo 

Substituting  with  equation  (28)  and  evaluating  the  resulting  integrals  using 
equation  (A-1),  we  find 


<^n  —  +  03)  • 

Similarly,  it  can  be  shown  that 


*^22  ~  (®1  +  +  03)  ) 

0-33  =  ^  («1  +  “2  +  803)  . 


(34) 


(35) 

(36) 


4.3  ID  Spectral  Densities 


First  let  us  consider  the  case  in  which  f  =  j  =  n  in  equation  (12).  Setting 
k2  =  k±  cos  9  and  ks  =  k±  sin  9  in  equation  (33),  we  find 


4>ii  (k)  =  Ai  {k)  +  A2  {k) 


kikj_cos"‘9  kikj_sin‘^  9 


k^ 


+  As  (k) 


(37) 


After  converting  equation  (12)  to  cylindrical  coordinates  and  performing 
the  integration  over  9,  we  have 


©11,1  {ki)  =2tt  [ 
Jo 


k\Ai  {k)  +  -klA2  {k)  +  -fcf  A3  {k) 


^dk^.  (38) 
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When  Ai  (k)  is  given  by  equation  (28),  the  terms  in  equation  (38)  can  be 
easily  integrated  using  equation  (A-2).  The  result  is 


FL(k)  = 


B (1/2, 1/3) 


4ai£i 


a2k^£2 


a3k^£n 


5  (1  +  A:2^2)5/6  6(1 +  A:2£2)  11/6  6(1  + A:2£2)  11/6 


(39) 


Next  let  us  find  an  expression  for  ©22,1  (k)-  First  note,  from  equation  (23), 
that  ‘F22  (k)  can  be  written 


kfkj^cos^O  kf  +  2kfkj^sm^  0  +  kj 

$22  (k)  =  (k)  - +  yl2  (k)  ^ - 1 

.  ^  /cl  (cos^  0  —  cos^  0) 

+  ^3  (k)  - ri - 


sin*^  0 


(40) 


Integrating  with  respect  to  0,  we  find 


©22,1  (ki)  =  2tt  [ 
Jo 


^kjk\Ai  (k)  +  +  klk'i  +  ^/ci^  A2  (k)  +  (k) 


(41) 


The  result  of  the  integration  over  k_\_  is 


FT(k)  = 


1 


aik'^£\ 


5(1/2,  2/3) 
02^2^1 


+ 


302/^2 


+ 


6(1 +  A:2£2)  11/6  10(1 +  A:2£2)  5/6 

lla2k^£2  ,  as/’s 


+ 


+ 


3(1  +  A:2£2)11/6  18(1  +  A:2^2)17/6  IQ  (1  +  A;2£2)5/6 

Similarly, 


(42) 


Fv(k)  = 


1 


aik‘^£\ 


5 (1/2, 2/3) 
303^3 


+ 


02/^2 


+ 


6(1  +  A:2£2)11/6  io(1  +  A:2£2)5/6 

a3k‘^£\  lla3k‘^£^ 


+ 


+ 


10(1  +  A:2£2)  5/6  3(1  +  fe2£2)  11/6  18(1  +  A:2£2)  17/6 


(43) 


In  the  inertial  subrange  {k£i  1),  it  is  widely  known  that  the  ID  spec¬ 
tra  satisfy  isotropy  and  are  proportional  to  A:“5/3  (Xhis  was  originally  sug¬ 
gested  by  Kolmogorov  (1941).)  In  the  inertial  subrange,  equations  (39)  and 
(42)  reduce  to 
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011,1  (A:)  ^ 


1 


5(1/2, 1/3) 
17 


155(1/2,1/3) 


^2^-2/3^-5/3^ 


and 


©22,1  (k)  = 


- - -  (- 

5  (1/2, 1/3)  Ve 

68 


455(1/2,1/3) 


+  l„3f3-/=)  fc-V3 

^2^-2/3^-5/3_ 


(Eq  (32)  was  used  in  deriving  the  equations  above.)  Note  that  both 
spectra  have  the  proper  dependence.  Furthermore,  the  ratio 

©11,1  {k)  /©22,i  (fc)  equals  4/3,  another  well-known  result  for  the  inertial 
subrange  (Panofsky  and  Dutton,  1984). 

The  ID  spectra  for  the  top-down  von  Karman  model  are  compared  to 
Kaimal  et  al's  equations  in  figures  3  to  5.  In  making  these  plots,  I  selected  the 
parameters  in  the  manner  to  be  discussed  in  section  4.5.  Agreement  is  fairly 
good  for  the  stream  wise  spectrum  F^.  The  peak  in  the  von  Karman  model 
is  somewhat  sharper,  and  the  spectral  levels  for  small  kz  are  lower.  Similar 
comments  apply  to  the  cross-stream  spectrum.  Ft,  with  the  mismatch  for 
small  kz  being  even  more  severe.  The  vertical  spectrum  Fy  is  rather  pe¬ 
culiar.  It  has  two  peaks,  one  at  kz  ~  0.3,  and  the  second  at  kz  ~  10.  This 
behavior  is  quite  unphysical  and  does  not  agree  well  with  the  Kaimal  et  al's 
spectra.  The  behavior  is  an  unfortunate  difficulty  with  the  top-down  ap¬ 
proach  to  modeling  the  anisotropic  spectra:  it  can  result  in  multiple-peaked 
ID  spectra. 


4.4  2D  Correlation  Functions 

From  equation  (23),  we  see  that 

^ii(0,A:2,fe3)  =  ^i(A:±),  (44) 

where  kj_  =  k2  +  k^.  Hence  equation  (14)  reduces  to 

bu,i  (r±)  =  /  (^-l)  •  >^-l)  (45) 

Rewriting  the  integral  in  cylindrical  coordinates  and  performing  the  inte¬ 
gration  over  angle,  we  have 

roo 

bu,i{r±)  =  27r  Ai  {k±)  Jq  {k±r±)  k±  dk±,  (46) 

40 
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where  {■)  is  the  Bessel  function  of  the  first  kind. 

Now,  substituting  equation  (28)  into  equation  (46),  we  have 


bn,i  {r±) 


3ai 

25(7/2,1/3) 


klel 


Jo  {k±r±)  dk±. 


(47) 


To  perform  the  integrations,  we  can  use  the  recurrence  formula  for  Bessel 
functions 

OiTl 

Jn+l  —  Jn  (x)  Jn—l  (x)  ■ 

X 


Hence 


k\Jo  {kxrx) 


^  J2  (fc±r_L)  -  J3  (/c_Lr_L)  +  fei  J4  {kxrx)  ■ 
r  I  r  I 


^11,1  (?'±) 


Now  the  integration  over  k±_  can  be  accomplished  using  equation  (A-3). 
The  result  is 


4aif  1  (  r  8^1  /r_L  \ 

5V^r(l/3)  V24/  [ri 


_S^K- 

r_L 


1/6 


+  -^7/6 


(48) 


where  Ky  (•)  is  the  modified  Bessel  function  of  the  second  kind.  The  for¬ 
mula  for  622,2  (?’±)  is  the  same,  except  with  02  replacing  ai,  and  ^2  replacing 

h. 

To  find  the  approximate  form  of  the  2D  correlation  for  small  separations 
{r±/l  <C  1),  we  use  the  following  series  expansion  for  the  Bessel  functions: 


m 


r(i-v)/cr  r(z/) 

2v  V2y  2(1-1/) 


The  result  of  substituting  this  series  into  equation  (48)  is 


&ii,i  (J’x)  - 


4r  (5/6)  ai^i 

50fr  (1/3) 


187r(l/6)  /^\5/3  2- 

6or(5/6)  \2eJ  ^  \2eJ 


Hence  the  2D  structure  function,  equation  (15),  is 


dii,i  {r±)  ~ 


lir(l/6)  2  5/3 

30V^r(2/3)  “  ^ 


36r(5/6)ai£i 
5^r(i/3)  ’ 


(49) 


(50) 


(51) 


where 


34cj2r(2/3) 

5r(i/3) 


{2i)~^^^ 


(52) 


is  the  structure-function  parameter.  The  first  term  in  equation  (51)  is  the  one 
for  isotropic,  inertial-subrange  turbulence.  It  is  the  most  important  term  for 
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very  small  r±.  But  as  r±  approaches  i  and  ii,  the  other  terms  quickly  be¬ 
come  important.  Hence  inertial-subrange  approximations  for  the  2D  struc¬ 
ture  function  generally  are  unsatisfactory. 

Note  that  the  2D  structure  functions  for  the  top-down  anisotropic  von 
Karman  model  are  independent  of  the  direction  of  the  displacement.  This 
is,  in  fact,  a  consequence  of  using  the  simplified  form  of  anisotropy  im¬ 
plied  by  equation  (23).  It  is  not  a  general  feature  of  anisotropic  turbulence. 
On  the  other  hand,  the  model  predictions  do  depend  on  the  orientation  of 
the  velocity  fluctuations.  For  example,  dny  {r±)  /  ^22, 2  i'f'x)- 

The  2D  structure  function  for  the  along-wind  direction  (dii,i  (rj,))  for  the 
top-down  anisotropic  von  Karman  model  is  plotted  in  figure  9.  The  corre¬ 
sponding  structure  function  for  crosswind  direction  (^22, 2  {r±))  is  plotted 
in  figure  10.  These  predictions  were  made  using  the  parameter  selections 
discussed  in  section  4.5.  The  model  predicts  that  the  2D  structure  function 
is  larger  in  along-wind  direction  than  in  the  crosswind  direction.  This  is  to 
be  expected,  since  the  variance  of  the  along-wind  velocity  fluctuations  is 
higher  than  for  the  crosswind  fluctuations. 

4.5  Parameter  Selection 

In  this  section,  I  determine  the  parameters  by  forcing  the  anisotropic  von 
Karman  model  to  have  the  same  variances  and  TKE  dissipation  rates  as  the 
empirical  equations  developed  by  Kaimal  et  al  (1972).  First  let  us  consider 


Figure  9.  Predictions  of  2D  structure  fimction  for  along-wind  direction,  for  Kologorov 
model,  various  von  Karman  models,  and  Kristensen  et  al  model. 
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Figure  10.  Predictions  of  2D  structure  function  for  crosswind  direction,  for  Kolmogorov 
model,  various  von  Karman  models,  and  Kristensen  et  al  model. 


the  variances.  By  comparing  equations  (34)  to  (36)  to  table  1,  we  have  the 
following  matrix  equation: 


’o.8  0.1  O.r 

ai 

’4.77' 

0.1  0.8  0.1 

02 

= 

2.68 

.0.1  0.1  0.8. 

.  03  . 

.  1.46. 

Solving  this  system  yields  the  values  for  the  a/s  shown  in  table  2. 

The  parameters  li  can  be  determined  by  matching  the  TKE  dissipation  rate 
in  the  inertial  subrange.  For  Mi  ^  1,  equation  (39)  becomes 


011,1  (A:)  ^ 


17 


155(1/2,1/3) 


-2/3, -5/3 


But  Kaimal  et  al  (1972)  have  demonstrated  the  following  equation 


(53) 


where  e  is  the  TKE  dissipation  rate,  and  oi  is  a  constant  approximately 
equal  to  0.56  (Kaimal  et  al,  1972,  Kristensen  et  al,  1989).  Therefore 


17 

155(1/2,1/3) 


^2^-2/3 


(/I  2/3 
2 


(54) 
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Table  2.  Values  for  top-down  von  Karman  model  parameters  Ui  and  £i. 


i 

ai 

li 

1 

5.54m^ 

4.922 

2 

2.56m^ 

1.542 

3 

0.813m* 

0.2772 

It  can  be  shown  that  the  TKE  dissipation  rate  in  a  neutral  shear  layer  is 

e  =  (55) 

KZ 

where  u*  is  called  the  friction  velocity,  z  is  the  height  from  the  ground,  and 
K  ~  0.4  is  called  von  Kdrmdn's  constant  (Stull,  1988).  Solving  equation  (54) 
for  i,  we  have 


15aiB(l/2,l/3) 


KZ  „  a/  ; 

^  ~  0.378^ 
ui  ui 


~  4.92z. 


The  parameters  £2  and  £3  now  follow  easily  from  equation  (29).  The  result¬ 
ing  values  for  the  £i  are  shown  in  table  2. 
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5.  Von  Karman-Based  Anisotropic  Model,  Bottom-Up  Approach 


Since  our  initial  attempt  to  derive  a  simple  anisotropic  von  Karman  model 
was  not  very  successful  in  producing  ID  spectra  that  agree  well  with  data 
(sect.  4),  we  try  in  this  section  a  second  approach  similar  to  the  one  devised 
by  Kristensen  et  al  (1989).  The  idea  is  to  use  the  same  set  of  equations  for  the 
ID  spectra  as  in  the  isotropic  von  Karman  model  (although  with  different 
parameters),  so  that  the  isotropic  and  anisotropic  versions  match  when  the 
variances  and  length  scales  are  all  equal. 
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equation  (27)  in  the  Wilson  report,  which  gives  the  transverse,  ID  spectrum 
in  the  isotropic  von  Karman  model.  If  we  define 


^  (i .  rrw) 

we  can  conveniently  state  equations  (57)  to  (60)  as 

Fi(/c)  =  Fo  (fill,  £i,0, 1/3; /c) , 

FT{k)  =  ^Fo(cj22,^2,5/3,1/3;/c), 

and 

Pv  (k)  =  ^Fo  (0-33,  £3, 5/3, 1/3;  k)  . 


(61) 


(62) 

(63) 

(64) 


The  ID  spectra  for  the  bottom-up  von  Karman  model  are  plotted  in  figures 
3  to  5.  Agreement  with  the  Kaimal  et  al  equations  is  fair.  The  main  difference 
for  all  three  cases  is  that  the  von  Karman  model  curves  are  more  sharply 
peaked  than  the  Kaimal  et  al  equations  and  have  less  energy  at  small  kz. 


5.2  Anisotropic  Energy  Spectra 


To  derive  the  anisotropic  energy  spectra,  we  follow  the  method  developed 
by  Kristensen  et  al  (1989).  Their  equations  (49)  to  (50)  for  the  A*  {k)  can  be 
written  in  the  following  form: 


Al  {k) 
A2  {k) 


k  d  1  dFi  {k)  1 

dvr  dk  k  dk  27r 

k  d  1  dFi  {k)  1 

Att  dk  k  dk  27r 


-/3 

a  —  3/3  —  C 


(65) 

(66) 


and 


7I3  {k) 


k  d  1  dFi  {k)  1 
Att  dk  k  dk  27r 


—a  (ji  —  3/3  (jt 


(67) 


The  first  term  in  equations  (65)  to  (67)  represents  the  isotropic  contribution. 
By  differentiation  of  equation  (62),  this  term  is  found  to  be 


k  d  IdFL  jk)  _  55cjfi _ k^if 

Att  dk  k  dk  367ri3  (1/2, 1/3)  (^  _l_ 


(68) 


Derivation  of  the  terms  in  equations  (65)  to  (67)  that  involve  the  functions 
a  (s),  /3  (s),  and  (  (s)  begins  by  defining  the  "residual"  functions 


H  (k)  =  Ft  {k)  -  Fv  {k) , 


(69) 
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and 


J  {k)  =  Fl  (k)  -  kF'L  {k)  -  [Ft  (k)  +  Fv  (k)] .  (70) 

In  isotropic  turbulence,  the  residual  functions  are  zero.  It  is  convenient  also 
to  define 

f  (s)  =  H  (71) 

and 

g{s)  =  j[s-^/^)s^  (72) 

where 

s  =  k-^.  (73) 

From  equations  (62)  to  (64),  we  find 

2 

/(s)  =  y  [Fo(a22,^2,5/3,l/3;s-i/2)  -Fo(u33,4,5/3,1/3;s-i/2^]  ,  (74) 

and 

9{s)  =  y  2Fo ^1,5/3, 1/3;  -Fo  1^0-22, ^2, 5/3, 1/3; 


-  Fo((T33,4,5/3,1/3;s-^/2^  .  (75) 

The  functions  a  (s),  (3  (s),  and  C  (s)  in  equations  (65)  to  (67)  can  be  deter¬ 
mined  from  the  integrals  (Kristensen  et  al,  1989) 

«  (s)  =  ^  (t)  dt  -  — ^  ^  ^l+l/v^/"  (i)  dt,  (76) 

„l/3  ps 

[3  (s)  =  —  /  (t)  dt,  and  (77) 

3  Jo 

C  (s)  =  -s~^  f  fg'"  (t)  dt.  (78) 

Jo 

Hence  we  need  to  solve  integrals  of  the  form 


I  {a,i,q;s)  =  ^-  t^  f^Fo  (a,  5/3, 1/3;  f  dt,  (79) 

and  we  have 
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a(s)  =  I  (^a22,i2,-^;s^  -  I  (^a33,e3,^-,s^  +  I  (^a22,h,--^;s^ 


-  I  (  (T33,£3, S 


P{s)  =  21  (^an,h,^-,s^  -  I  (^a22,i2,^;s^  -  I  (^a33,£3,^;s^ 


(80) 


(81) 


and 


C(s)  =  21  {an,h,-l;s)  -  /  ((T22,  ^2, -1;  s)  -  /  (0-33,  £3, -1;  s)  .  (82) 

It  can  be  shown  that 


f  „  -i/2\  v  (i' —  a  +  1/2)  1“^  ( 

{a,e, 


a  (1/  +  3/2) 
u-a  +  l/2' 


from  which  follows 

r 


dt^  L 


t^Fo{a,£,a,v-t  =  C3^Fo  (a,  03, 2/3;  f  , 


where 


“1“  1/2)  Q-n  {^n  F  8/2) 

Cn+1  —  ,  -I  /r,  Cn,  Q-n+1  —  ;  777T ;  ^n+1  —  +  1; 


+  1/2 


^'n  ~  Q-n  +  1/2 


and 


Co  —  1,  oo  —  «  —  g,  t'o  —  —  2' 


(The  Cj's  defined  in  this  section  are  different  from  those  in  sect.  3  .)  It  can  be 
shown  that 

140  23  10 


C3  27  ’  “3  ^5  >  '^s  2 


Hence 


^  r. 

dF  - 


f2Fo(a,A5/3,l/3;ri/2) 


27  ’  15’  3 


4675 


(t/f) 


-1/6 


216  (1/2,  1/3)  (1  +  / /£2)23/6  15  1+  t/(2 


1  - 


23  1 
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With  the  help  of  equation  (A-7),  we  now  find 


I{a,i,q]s)  = 


4675  /sy 

'4^5(l/2,l/3)  \j^) 

-B 


23  fll 

-  9, 3  +  <? 


:/£2/(i+s/<;2)  (^y  -  ^,2  +  g 


(83) 


where 

B,{a,b)=  r  s^^-^il-sf-^ds  (84) 

Jo 

is  the  incomplete  beta  function.  Substitution  of  this  result  into  equations  (80) 
to  (82)  yields  the  functions  a{s),  (3  (s),  and  C{s),  and  hence  the  anisotropic 
energy  spectra. 

The  anisotropic  energy  spectra  for  the  bottom-up  von  Karman  model  are 
shown  in  figures  6  to  8.  One  notable  feature  of  the  As-curve  is  that  it  is  ini¬ 
tially  negative  and  then  becomes  positive  for  kz  >  2.  This  behavior  is  some¬ 
what  alarming,  since  an  energy  spectrum  must  strictly  be  non-negative. 
The  behavior  could  result  from  vertical  inhomogeneities  characteristic  of  a 
shear  layer,  which  are  neglected  in  the  present  model. 


5.3  2D  Correlation  Functions 


The  2D  correlation  function  6ii  i  (r^)  can  be  found  by  substituting  equation 
(65)  into  equation  (46).  The  isotropic  contribution,  given  by  equation  (68), 
is  therefore 


roo 

'k^  d  1  dFL{k)' 

Jo  {kr) 

dk  - 

POO 

k^k\ 

L 

2  dk  k  dk 

185(1/2,1/3), 

Jo 

.(1  +  A:2£2)17/6_ 

Jo  (kr)  dk. 


(85) 


The  recurrence  relation  for  Bessel  functions  allows  us  to  write  Jq  {kr)  = 
{2/kr)Ji  {kr)  —  J2  {kr).  This  leaves  two  integrals,  each  of  which  can  be 
solved  using  equation  (6.565.4)  in  Gradshteyn  and  Ryzhik  (1994).  The  re¬ 
sult  is 


k^  d  1  dPL  {k) 
2  dk  k  dk 


Jo {kr)  dk 


0fr(i/3) 


r 

-  -^K 


1/6 


(86) 


Calculation  of  the  anisotropic  contribution  to  6114  (r)  is  rather  difficult.  We 
need  to  solve  integrals  of  the  form 

I  {a,  i,q]r)  =  j  I  (^a,  i,  q;  do  {kr)  k  dk,  (87) 
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in  which  I  (a,  i,  q;  k~‘^)  contains  incomplete  beta  functions  (eq  (83)).  One 
possible  approach,  based  on  rewriting  the  incomplete  beta  function  as  a 
Meijer's  G-function,  is  discussed  in  the  appendix.  Using  equation  (A-12), 
we  find 


I  {a,£,q;r) 


qa'^i  f  ^  \ 

2V^r(i/3)  \Yi) 


17 jQ,  2,  23/6  -  q,  17/6 
17/6 -g,  35/6 


17/6,  2,  23/6 -g,  17/6 
17/6  -  q,  29/6 


(88) 


^11,1  (r) 


For  the  2D  correlation  function  in  the  along-wind  direction,  we  now  have 
from  equation  (65) 


2(Tn£i 

V^r(i/3) 


r 

- K 

2ii 


1/6 


Pir)-C  (r) , 


(89) 


where,  from  eq  (27), 

]5{r)  =  2l  -7  (^(T22,^2,  -7  (^0-33,4,  •  (90) 

Equations  for  622,2  (r),  633^3  (r),  a  (r),  and  ( (r)  follow  with  obvious 
replacements. 

Since  routines  for  computing  Meijer's  G-functions  are  not  widely  avail¬ 
able,*  calculating  the  2D  correlations  using  equation  (88)  is  not  normally 
practical.  Rather,  they  must  be  calculated  by  performing  the  integration 
in  equation  (87)  numerically.  Generalization  of  the  isotropic  von  Karman 
model  to  the  anisotropic  case,  using  a  bottom-up  approach,  therefore  pre¬ 
vents  us  from  obtaining  important  results  such  as  the  2D  correlations  in  a 
convenient  analytical  form — a  rather  disappointing  outcome. 

Computations  of  the  2D  structure  functions,  found  by  numerical  integra¬ 
tion,  are  shown  in  figures  9  to  10.  The  results  are  observed  to  differ  signifi¬ 
cantly  from  the  top-down  approach. 


5.4  Parameter  Selection 

As  was  the  case  for  the  top-down  approach,  I  choose  the  variances  for  the 
bottom-up  anisotropic  von  Karman  model  to  match  Kaimal  et  al's  results, 
shown  in  table  1.  The  length  scale  parameters  are  chosen  to  produce  the 

*H.  Auvermann  has  shown  that  the  commercial  software  package  Mathematica  contains 
a  routine  to  compute  Meijer's  G-functions. 
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correct  inertial-subrange  asymptote,  just  as  in  section  4.5.  From  equations 
(53)  and  (108  ),  we  have  for  Mi  »  1, 


Fl  {k)  ^ 


B (1/2, 1/3) 


Solving  for  ii,  we  find 


= 


20-2 


11 


aS  (1/2, 1/3) 


-|3/2  3 

^  c.  0.350^^  c.  3.64Z. 


It  follows  from  equation  (58)  that 


Ft  (k)  =  Fv  (k)  =  -Fl  (k) 

in  the  inertial  subrange.  Hence  for  ki2  ^  1 , 


and 


Similarly, 


i2  = 


2  <722 


-,3/2 


aS  (1/2, 1/3) 


KZ 


~  0.350- 


.0-22^ 


~  1.53z. 


4  = 


^*^33 

aS  (1/2, 1/3) 


KZ  O-ooZ 

^  ~  0.350^1-  ~  0.617Z. 


ui 


ui 


(91) 


(92) 


(93) 


(94) 
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6.  Kristensen  et  al  Model 


Krasnenko  et  al  (1997)  have  previously  suggested  using  the  Kristensen  et  al 
(1989)  model  to  calculate  statistics  of  acoustic  field  fluctuations.  This  model 
is  very  similar  to  the  bottom-up  von  Karman  model  discussed  in  the  pre¬ 
vious  section.  In  fact,  I  based  the  derivation  in  the  previous  section  largely 
on  Kristensen  et  al's  model.  The  main  comparative  advantage  of  the  Kris¬ 
tensen  et  al  model  is  its  flexibility.  It  contains  additional  parameters  that 
allow  one  to  tailor  the  shape  of  the  spectral  "peak"  to  the  data  set.  The 
disadvantages  are  that  the  equations  are  somewhat  more  complicated  and 
that,  in  some  cases,  deriving  values  for  all  of  the  parameters  may  be  diffi¬ 
cult.  The  equations  for  the  Kristensen  et  al  model  are  summarized  below  in 
a  slightly  different  form  than  in  the  original  paper.  The  changes  hopefully 
make  computations  somewhat  simpler. 


6.1  ID  Spectral  Densities 

Kristensen  et  al  used  the  following  equations  for  the  ID  spectra  for  wave 
numbers  corresponding  to  the  direction  of  the  mean  wind: 


FL{k) 


Frik) 


011,1  (A:)  = 


©22,1  (k)  = 


Lii^crh 


TT 


1  + 


fcLi 


L  o.{ai) 


2^1  f  ’ 


L22, 1(^22 


^  3 


kL22: 


a{ti2) 


2fM2 


‘l-K 


1  + 


kL22,l 


a{ti2) 


(95) 


(96) 


and 


Fv  {k)  =  033,1  {k) 


-^33,10'33 

2tt 


fc^33,l 

“(ms) 


2M3 


fc^33,l 

“(ms) 


|^5/6/i3+l  ’ 


(97) 


where 


a{n) 


TT/r 


(98) 


The  parameters  m  control  the  sharpness  of  the  spectral  peak;  larger  values 
give  a  sharper  peak.  A  value  of  1  corresponds  to  the  bottom-up  anisotropic 
von  Karman  model. 
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The  ID  spectra  for  the  Kristensen  et  al  model  are  plotted  in  figures  3  to 
5  and  compared  to  Kaimal  et  al  equations,  as  well  as  the  two  versions  of 
the  von  Karman  model.  (Selection  of  the  parameters  for  this  model  will  be 
discussed  in  sect.  6.4.)  It  is  seen  that  the  Kristensen  et  al  model  has  enough 
flexibility  in  its  equations  to  reproduce  Kaimal  et  al's  empirical  equations 
quite  well. 


6.2  Anisotropic  Energy  Spectra 

Derivation  of  the  anisotropic  energy  spectra  A*  (k)  proceeds  much  as  it  did 
in  section  5,  except  that  the  equations  become  somewhat  more  complicated 
because  of  the  more  complicated  form  for  the  ID  spectra.  The  first  term  in 
equations  (65)  to  (67),  representing  the  isotropic  spectra,  is  found  by  differ¬ 
entiation  of  equation  (96)  to  be 


k  d  IdFiik) 
47r  dk  k  dk 


5cjfi 


Pmi-2^2m2+i 


(pi) 


1  + 


.  o-iai)  \  j 


5/6A41+2 


2pi  -|- 


11  kLiii 

3  [a(/ii) 


(99) 


The  remaining  terms  in  equations  (65)  to  (67),  involving  the  functions  a  (s), 
(3  (s),  and  ( (s),  represent  the  anisotropic  contributions.  Since  the  reader  can 
refer  to  the  solution  for  these  functions  in  Kristensen  et  al  (1989),  it  will  not 
be  repeated  here.  The  functions  are  given  by  the  equations 


“(■s)  =  70  0-22,^22,1,  -70  ^^,0-33, -^33, 1,^3; 

+  70  (“^’'^22,  -^22,1,  M2;  -  70  (“^’^33,  .^33,1,  M3;  , 


(100) 


P{s)  —  270  (Til,  L114,  pi;  -70  ^-,0-22,^22,1,  M2;  -70  ^2,0-33,  .b334,p3;s^  ,  (101) 


and 


C  (s)  —  270  (  — 1, 0-11,  L117,  Mi;  s)  —  7o  (  —  1, 0-22,  ^22,1,  M2;  s)  —  7o  (  —  1,0-33,  1^33,1,  M3;  s)  .  (102) 

The  function  70  is  in  turn  given  by 

^  ^  W  Q  _  1  \ 

70 {9, tt)  • 
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where 


(104) 


and 


T]  = 


a  ifi) 
£2 


Cl  ifi)  =  40  (1  -  fi)  (1  -  2//)  (2  -  fi) 

C2  (m)  =  ^  (1  -  m)  (1  -  2li)  (5  +  6/i) 

C3  (m)  =  f  (5  +  6m)  (5  +  12m)  (7  -  12m) 
C4  (m)  =  i  (5  +  6m)  (5  +  12m)  (5  +  18m)  . 


The  anisotropic  energy  spectra  for  the  Kristensen  et  al  model  are  plotted  in 
figures  6  to  8  and  compared  to  the  models  discussed  in  previous  sections. 
The  curves  for  the  Kristensen  et  al  model  are  similar  to  the  bottom-up  von 
Karman  model,  although  they  appear  more  flattened.  As  was  the  case  for 
the  bottom-up  von  Karman  model,  the  Kristensen  et  al  model  exhibits  neg¬ 
ative  A3  for  small  kz.  Since  the  Kristensen  et  al  model  fits  Kaimal  et  al's  ID 
spectra  so  well,  the  cause  of  this  feature  must  either  be  that  the  ID  spectral 
measurements  are  inaccurate  or  that  the  turbulence  cannot  be  adequately 
modeled  using  equation  (23). 


6.3  2D  Correlation  Functions 

In  order  to  find  the  2D  correlations,  we  need  to  solve  the  integral  equation 
(46),  with  Al  given  by  equation  (65).  For  the  Kristensen  et  al  model,  we 
must  calculate  integrals  of  the  form 

70  (g,  tr,  L,  m;  Jo  (/cr) /c  dfc. 

I  have  been  unable  to  find  a  method  for  calculating  these  integrals.  Appar¬ 
ently  one  cannot  even  express  the  result  in  terms  of  Meijer's  G-functions, 
as  was  possible  for  the  bottom-up  von  Karman  model.  Fortunately,  finding 
the  2D  correlations  by  performing  a  numerical  integration  is  not  difficult. 

Predictions  from  the  Kristensen  et  al  model  are  compared  to  the  von  Karman 
family  in  figures  9  to  10.  Comparisons  to  the  anisotropic  Gaussian  and 
Mann  models  (to  be  discussed  later  in  this  report)  are  shown  in  figures  11 
to  12.  The  Kristensen  et  al  model  predicts  a  strong  enhancement  of  the  2D 
structure  function  in  the  along-wind  direction,  and  a  slight  enhancement  in 
the  crosswind  direction. 
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Figure  11.  Predictions  of  2D  structure  function  for  along-wrnd  direction  for  Kristensen  et  al, 
Gaussian,  isotropic  von  Karman,  Kolmogorov,  and  Mann  models. 


Figure  12.  Predictions  of  2D  structure  function  for  crosswind  direction,  for  Kristensen  et  al, 
Gaussian,  isotropic  von  Karman,  Kolmogorov,  and  Mann  models. 
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6.4  Parameter  Selection 


Kristensen  et  al  chose  the  variances  and  integral  length  scales  in  their  model 
to  match  those  found  by  Kaimal  et  al,  and  hence  the  values  for  these  param¬ 
eters  are  the  same  as  in  table  1.  They  also  determined  the  parameters  jji  by 
fitting  their  model  to  Kaimal  et  al's  equations.  The  result  is 

Pi  =  0.52  1 

P2  =  0.49  \  ■  (106) 

P3  =  0.68  j 

Note  that  the  p/s  deviate  significantly  from  1,  indicating  that  the  actual 
spectra  are  less  peaked  than  in  the  von  Karman  model. 
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7.  Gaussian  Model 


Although  it  is  rather  unrealistic,  the  isotropic  Gaussian  model  has  proven 
quite  popular.  The  main  advantage  of  the  isotropic  Gaussian  model  is  that 
the  equations  for  the  correlation  and  spectral  functions  are  easily  obtain¬ 
able.  In  this  section,  we  explore  how  the  Gaussian  model  can  be  gener¬ 
alized  to  the  anisotropic  case  using  a  bottom-up  approach.  Since  conve¬ 
nient  closed-form  results  could  not  be  obtained  for  quantities  such  as  the 
2D  correlations  in  the  bottom-up  von  Karman  model,  it  is  interesting  to  see 
whether  such  results  can  be  obtained  for  the  Gaussian  model. 


The  ID  spectra  for  the  anisotropic  Gaussian  model  are  plotted  in  figures  3 
to  5,  which  also  show  the  models  discussed  in  previous  sections  of  this 
report.  (Selection  of  the  parameters  for  the  Gaussian  model  will  be  de¬ 
scribed  in  sect.  7.4.)  The  Gaussian  model  is  observed  to  agree  with  Kaimal 
et  al's  equations  for  very  small  kz.  However,  the  spectral  peak  in  the  Gaus¬ 
sian  model  is  much  too  sharp,  and  the  inertial-subrange  decay  is  much  too 
rapid. 
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7.2 


Anisotropic  Energy  Spectra 

We  will  again  determine  the  anisotropic  energy  spectra  using  equations 
(65)  to  (67).  The  isotropic  term  in  those  equations  is,  for  the  Gaussian  model. 


k  d  1  dFi  (/c)  f 

dvr  dk  k  dk  327r^/2  ^  y  A  j 


(112) 


Using  the  same  method  as  in  section  5.2  to  calculate  a  (s),  (3  (s),  and  C  (s), 
we  arrive  at  the  residual  functions 


/{»)  =  y  L 


2 

^  To  (cr22,  ^2;  -  To  (0-33,  ^3;  ,  and 


o2  . 


5(^)  =  y  L 


2Ff 


—  -To  (^0-22,  ^2;  s  —  -To  (^(T33,b3;  S 


where 


Fo{a,e-,k)  = 


aH 

2A 


2ff2' 


1  + 


kH- 


2n2' 


exp 


kH 


(113) 


We  again  need  to  find  integrals  of  the  form  given  by  equation  (79),  although 
with  the  function  To  given  by  equation  (111),  we  can  show  that 


^  r 
dF  - 


t^Fo 


<jH 

2^ 


5^6 

+ 


64T  128T 


exp 


'At 


The  integrals  can  now  be  found  using  equation  (A-9),  with  result 


I  ((T,b,g;s)  = 


/4sV 


r(3  +  c?,^)  -|r  (2  +  g,^ 


4s 


4s 


(114) 


The  functions  a  (s),  (3  (s),  and  (  (s)  can  now  be  determined  using  equations 
(80)  to  (82). 

The  anisotropic  energy  spectra  for  the  Gaussian  model  are  plotted  in  fig¬ 
ures  6  to  8.  They  have  a  rather  strange  and  unphysical  appearance,  exhibit¬ 
ing  multiple  positive  and  negative  peaks  as  kz  is  varied. 


7.3  2D  Correlation  Functions 


It  is  not  difficult  to  calculate  the  isotropic  contribution  to  the  2D  correlation 
functions  in  the  Gaussian  model.  Substituting  equation  (112)  into  (46),  we 
have  the  following  integral  for  the  isotropic  contribution: 


k^  d  1  dFL  (fc) 
2  dk  k  dk 


Jo (kr)  dk 


Jo  (kr)  k^  dk. 


(115) 
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Using  the  recurrence  relation  for  Bessel  functions,  we  have  two  integrals, 
each  of  which  can  be  solved  using  Gradshteyn  and  Ryzhik's  (1994)  equa¬ 
tion  (6.63.4).  The  result  is 


d  1  dFi  (fc) 
2  dkk  dk 


Jo  {kr)  dk 


(116) 


The  calculations  for  the  anisotropic  contributions  are  more  complicated,  re¬ 
quiring  us  to  integrate  incomplete  gamma  functions.  Specifically,  we  must 
calculate  integrals  of  the  form 


/  {a,£,q;r) 


(a,e,q;k 


Jo  (kr)  kdk  = 


qa‘^1^  f°^  /  4  X”? 

Sy/ir  Jo  V  k‘^£‘^  ) 


T 


T 

2 


Jo  (kr)  k  dk. 


(117) 


To  solve  an  integral  of  this  form,  I  will  rewrite  the  incomplete  gamma 
functions  in  terms  of  confluent  hypergeometric  functions.  From  equa¬ 
tions  (13.1.33)  and  (13.6.28)  in  Abramowitz  and  Stegun  (1965),  we  have 


T  (a,  X)  =  e-"/2a;“/2-l/2^^/2-l/2,-a/2  {x)  ,  (118) 

where  FFk,/^  (•)  Whittaker's  function.  Hence 


I  {a,£,q;.s) 


/4sy/2  / 


X 


^l+g/2,3/2+g/2 


5 

2 


1/2 

W^l/2+g/2,l+g/2 


.  (119) 


The  integration  can  now  be  performed  using  Gradshteyn  and  Ryzhik's 
(1994)  equation  (7.672.2),  with  result 


I  {a,£,q;r)  = 


3qa^£ 


2^(1  -  q) 


2J2  I  4,1  -  q]  1,2  -  q] 


5  / 

~4  f  8, 1  “  O';  1)  2  —  g;  — ^ 


(120) 


where  2F2  (•)  is  a  generalized  hypergeometric  function. 

Although  we  have  managed  to  integrate  the  anisotropic  part,  the  result  is 
not  very  helpful  from  a  practical  standpoint.  Routines  are  not  commonly 
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available  to  calculate  generalized  hypergeometric  functions.  Hence  the  2D 
correlations  must  be  calculated  numerically,  by  integrating  equation  (46). 
This  is  a  rather  surprising  outcome:  even  the  simple  Gaussian  model,  when 
generalized  to  the  anisotropic  case,  eludes  a  convenient  solution  for  the  2D 
correlation  function. 

Predictions  of  the  2D  structure  function  from  the  anisotropic  Gaussian 
model  (calculated  by  numerical  integration)  are  compared  to  the  Kris- 
tensen  et  al  model  in  figures  11  to  12.  The  Gaussian  model  converges  on  the 
Kristensen  et  al  model  for  sensor  separations  much  larger  than  the  height, 
since  the  two  models  are  essentially  the  same  at  large  scales  of  the  energy 
subrange.  However,  at  smaller  separations,  corresponding  to  the  inertial 
subrange,  the  Gaussian  model  predictions  are  quite  unrealistic.  Further¬ 
more,  the  Gaussian  model  exhibits  unphysical  peaks  in  the  2D  structure 
function  around  the  transition  between  the  inertial  and  energy  subranges. 


7.4  Parameter  Selection 


Since  the  Gaussian  model  does  not  have  a  realistic  inertial  subrange,  it  is 
logical  to  select  its  parameters  to  obtain  the  best  possible  fit  to  spectral  char¬ 
acteristics  of  the  energy  subrange,  such  as  variances  and  integral  length 
scales.  Hence  I  take  the  approach  here  of  selecting  the  Gaussian  model  pa¬ 
rameters  to  reproduce  the  variances  and  integral  length  scales  given  in  ta¬ 
ble  1.  To  reproduce  the  variances,  of  course,  we  simply  assign  the  variances 
in  the  Gaussian  model  to  those  given  in  the  table.  The  length  scales  follow 
from  equation  (18).  Evaluating  equations  (108),  (111),  and  (111)  at  /c  =  0,  we 
find 


-^22,1  —  -^33,1 


The  resulting  Gaussian  model  parameters  are  shown  in  table  3. 


Table  3.  Values  for  variances  and  length  scales  ii  corresponding  to  Gaussian  model. 


% 

2 

ii 

1 

4:.77ul 

6.212 

2 

2.68ui 

3.592 

3 

lA6ul 

0.7742 
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8.  Mann  Model 


8.1  Rapid  Distortion  Theory  and  Eddy  Lifetimes 

The  Mann  model  is  the  result  of  a  "first-principle"  theoretical  treatment  for 
turbulence  in  a  constant  shear  layer.  By  constant  shear,  it  is  meant  that  the 
gradient,  dU / dz  (where  U  is  the  mean  wind  speed),  is  constant.  Constant 
dU/dz  implies,  of  course,  that  U  is  proportional  to  z.  In  an  actual  shear 
layer,  U  is  proportional  to  the  logarithm  of  z.  Except  for  this  shortcoming, 
the  Mann  model  is  probably  quite  realistic. 

The  main  hypothesis  upon  which  Mann's  model  is  based  is  called  rapid  dis¬ 
tortion  theory.  The  idea  is  basically  this:  In  a  turbulent  shear  layer,  new  ed¬ 
dies  are  perpetually  being  created,  subjected  to  the  forces  of  shear,  and  then 
eventually  broken  down.  Between  the  times  when  an  eddy  is  "created"  and 
when  it  is  "destroyed"  by  the  shear,  there  is  a  time  at  which  the  eddy  exists 
in  a  state  that  is  "typical"  for  eddies  having  its  same  spatial  dimensions. 
Let  us  call  this  "typical"  time  the  eddy  lifetime.  This  argument  suggests  that 
we  can  derive  a  spectral  model  by  introducing  turbulence  with  some  en¬ 
ergy  spectrum  E  (k)  at  time  t  =  0  and  then  subjecting  the  turbulence  to 
shear  (i.e.,  rapidly  distorting  it)  for  a  time  interval  equal  to  the  eddy  life¬ 
time.  Then  the  spectrum  should  be  roughly  the  same  as  turbulence  in  the 
actual  shear  layer. 

The  initial  energy  spectrum  E  (k)  in  Mann's  model  is  the  von  Karman  en¬ 
ergy  spectrum,  equation  (1).  The  eddy  lifetime  being  given  in  nondimen- 
sional  form  (depends)  on  the  scale  of  the  eddies  by* 


r{ki) 


-,—;%-{k£)-^ 
3’  6  ’  3’  ^  ^ 


(121) 


in  which  T  is  a  model  parameter.  Values  for  T  in  the  range  of  2.6  to  3.8 
appear  to  be  reasonable  for  atmospheric  turbulence  (Mann,  1994). 

A  full  discussion  of  rapid  distortion  theory  and  Mann's  model  is  beyond 
the  scope  of  this  report.  The  results  of  the  analysis,  as  they  pertain  to  acous¬ 
tical  propagation  modeling,  will  be  summarized  in  the  following  sections. 

*The  nondimensional  eddy  lifetime  is  dimensionalized  by  multiplication  by  {dU/dz)~^. 
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8.2  3D  Spectral  Densities 

Mann's  equations  for  the  spectra  in  a  uniform  (constant  gradient)  shear 
layer  are 


$11  (k)  = 

E{ko) 

47r/cQ 

ko  —  k\  —  2/ciA:3oCi  +  {ki  +  k^  > 

(122) 

4^22  (k)  = 

E{h) 

dvr/cQ 

kl-kl-  2/c2A;3oC2  +  {kl  +  Cl  > 

(123) 

$12  (k)  = 

E{ko) 

47r/cQ 

—kik2  —  A:iA:3oC2  —  k2kzoCi  +  {kl  +  k^  C1C2  , 

(124) 

$13  (k)  = 

E{ko) 

Airk^k"^ 

—kik^Q  +  (^1  +  k^  Cl  ) 

(125) 

4^23  (k)  = 

E{ko) 

Airk^k"^ 

-k2k3o  +  {kl  +  C2  ,  and 

(126) 

4^33  (k)  = 

E{ko)  , 
Ank^  ' 

[kl  +  kl)  . 

(127) 

The  initial  (before  the  onset  of  shear  distortion)  wave  number  is 

ko  =  {ki,k2,k3o) , 

(128) 

where 

kso  =  k3-  fdki. 

(129) 

The  quantities  Q  are  given  by  the  equations 

Ci  =  Ci-^C2,  C2  =  ^Ci  +  C'2, 

kl  kl 

(130) 

where 

f3kl  {kl  -  2/c|o  +  l3kik3o) 

^  ~  k^  {kl  +  kl) 

(131) 

and 

k2kl  \(3ki{kl  +  klf^‘^ 

C2  = - 2——  arctan  — 2 . 

{kl  +  kl)^^  kl-Pk3oki 

(132) 

To  facilitate  implementation  of  Mann's  model  on  a  computer,  rewriting  the 
hypergeometric  geometric  function  in  equation  (121)  as  an  incomplete  beta 
function  is  helpful.  First,  by  applying  equation  (15.3.4)  in  Abramowitz  and 
Stegun  (1965),  we  have 


2F1 


'1  17_4_ 
.3’¥’ 3’ 
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Then,  by  applying  Abramowitz  and  Stegun's  equation  (6.6.8)  to  the  right- 
hand  side,  we  have 


2F1 


/I  5 


Hence  the  nondimensional  eddy  lifetime  can  be  expressed  as 


M 


B 


l/(l+fc2£2) 


1  5M-^/2 

3’  2 


(133) 


The  incomplete  beta  function  (with  argument  in  the  domain  [0, 1])  can  be 
computed  using  routines  available  with  software  packages  such  as  Matlab. 


8.3  Variances  and  Parameter  Selection 

The  variances  and  covariances  in  the  Mann  model  can  be  found  by  numer¬ 
ically  integrating  the  3D  spectra  (k),  in  accordance  with  equation  (9). 
The  results,  as  a  function  of  T,  are  shown  in  figure  13.  As  T  is  increased, 
the  variance  in  the  streamwise  component  increases.  So  too  does  the 
cross-stream  variance  0-22  increase,  although  at  a  lesser  rate.  The  vertical 
variance  0-33  decreases  slightly.  Note  also  that  the  model  predicts  a  negative 
covariance  af^,  between  the  streamwise  and  vertical  fluctuations,  which  is 
an  important  intrinsic  feature  of  a  shear  layer.  (All  of  the  models  discussed 
earlier  were  somewhat  unrealistic  in  that  they  predict  zero 

It  is  also  interesting  to  look  at  the  integral  length  scales  as  a  function  of  T. 
These  are  shown  in  figure  14.  (The  computation  involved  performing  the 
integration  in  eq  (12)  numerically  and  then  applying  eq  (18).)  When  T  = 
0,  we  have  the  usual  result  for  isotropic  turbulence  that  the  longitudinal 
integral  length  scales  (e.g.,  Tii,i)  are  twice  the  transverse  scales  (e.g.,  ^22,1 
and  T33  1).  All  of  the  integral  length  scales  shown  initially  increase  with 
increasing  T,  although  for  T  >  1.5,  T33  1  starts  to  decline. 

Maim  chose  the  model  parameters  using  a  least-squares  method  formu¬ 
lated  to  give  the  best  fit  with  ID  spectra.  In  this  report,  I  used  a  somewhat 
simpler  method  that  gives  similar  results.  The  main  idea  is  to  choose  T  to 
match  the  ratio  (y\\/(y\j,  measured  by  Kaimal  et  al.  From  table  1,  this  ratio 
should  be  4.77/1.46  =  3.27.  Using  the  data  shown  in  figure  13,  one  can 
show  the  value  of  T  that  produces  this  ratio  to  be  T  =  Tq  =  3.53.  Since 
Mann  (1994)  found  his  atmospheric  data  were  best  described  by  values  for 
T  in  the  range  of  2.5  to  3.5,  the  value  suggested  appears  reasonable. 
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Figure  13.  Variances  of  three  velocity  components  and  covariance  between  longitudinal  and 
vertical  velocities,  according  to  Mann  model.  (Co)variances  are  plotted  as  a  function  of  F,  a 
parameter  controlling  eddy  lifetime. 
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Figure  14.  Integral  length  scales  for  three  velocity  components  and  for  joint  fluctuations 
between  longitudinal  and  vertical  velocities,  according  to  Mann  model.  Length  scales  are 
plotted  as  a  function  of  F,  a  parameter  controlling  eddy  lifetime. 


(134) 


To  force  afi  equal  to  4.77u^  when  T  =  Tq,  we  need  to  set 

cj?  =  1.88u* 


in  the  initial  von  Karman  spectrum. 

In  order  to  reproduce  the  inertial-subrange  asymptote,  we  must  set 


£  = 


2r(5/6)ai 
y^r(l/3)ai 
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(135) 


8.4  2D  Correlation  Functions 


The  2D  correlation  functions  for  the  Mann  model  can  be  found  using  equa¬ 
tion  (14).  Since  (k^)  depends  on  the  orientation  of  k^,  the  double  in¬ 
tegral  cannot  be  reduced  to  a  single  one  involving  a  Bessel  function  as  in 
equation  (46).  However,  if  we  are  interested  in  sensor  separations  along 
just  one  of  the  orthogonal  coordinate  axes,  equation  (14)  can  be  simplified 
somewhat.  For  example. 


/OQ  poo 

/  4>ii  (0,/c2,A:3)  dfca 

-oo  —OO 


exp{ik2r2)  dk2,  (136) 


or 


/CXD  /*00 

/  4>ii  (0,/c2,A:3)  (iA:2 

-OO  — oo 


exp  (iksrs)  dks-  (137) 


The  advantage  of  these  equations  is  that  the  full  Fourier  transform  needs 
to  be  performed  along  just  one  axis.  A  simple  numerical  integration  can 
be  performed  along  the  other  axis.  Note  that  the  2D  structure  functions  in 
the  Mann  model  depend  on  the  direction  of  the  displacement.  This  distin¬ 
guishes  the  Mann  model  from  the  ones  considered  earlier. 

Predictions  of  the  2D  structure  function  for  the  Mann  model  are  shown 
in  figures  11  to  12.  Note  that  the  structure  functions  for  both  horizontal 
and  vertical  separations  are  nearly  the  same  for  small  separations  com¬ 
pared  to  the  height  (the  inertial  subrange),  diverge  for  moderate  separa¬ 
tions,  and  then  become  equal  for  large  separations.  The  reason  why  they 
become  equal  for  large  separations  is  that  bn,i  {r,  0)  and  (0,  r)  both 
approach  zero  for  large  r,  so  that  the  2D  structure  function  in  the  two 
cases  is  (r,  0)  ~  dn^  (0,r)  ~  26ii^i  (0,0).  For  large  separations,  the 
enhancement  of  the  along-wind  structure  function  dn^  predicted  by  the 
Maim  model  is  about  half  that  predicted  by  the  Kristensen  et  al  model.  The 
predictions  differ  even  more  in  the  crosswind  direction:  the  Mann  model 
predictions  for  d22,2  are  about  1/10  the  Kristensen  et  al  predictions. 
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9.  Anisotropy  and  the  Coherence  of  Sound  Waves 


So  far  in  this  report,  I  have  dealt  only  with  the  modeling  of  anisotropic 
turbulence,  without  discussing  its  significance  to  wave  propagation.  The 
effect  of  anisotropic  turbulence  on  the  coherence  of  a  propagating  sound 
wave  will  be  discussed  in  this  section. 


9.1  Mutual  Coherence  Function 


The  MCF  describes  the  coherence  between  a  pair  of  sensors  as  a  function 
of  their  separation  r,  normal  to  the  direction  of  propagation.  It  is  given  by 
(Wilson,  1997a) 

T  (r)  =  exp  — (r)  ,  (138) 

where  R  is  the  propagation  distance  (range);  k  =  27r//co  the  wavenumber, 
/  the  frequency,  cq  the  sound  speed;  and  du  (r)  the  2D  structure  function  of 
the  effective  index-of-refraction  fluctuations.  The  function  du  (r)  is  related 
to  its  counterpart  for  the  velocity  fluctuations  according  to  (Wilson,  1997a) 

du  (r±)  =  ^dii^i  (r_L)  ,  (139) 

Co 

where  the  direction  of  propagation  is  the  zth  axis,  and  is  a  displacement 
perpendicular  to  the  fth  axis. 

In  the  shear-layer  models  considered  in  this  report,  all  of  the  variances  scale 
in  proportion  to  u^,  and  all  of  the  length  scales  to  z.  From  equations  (13) 
and  (15),  we  see  that  da^i  (r^)  is  proportional  to  zu^.  Hence  let  us  define  the 
nondimensionalized  2D  velocity  structure  function  as 

=  (140) 

zui 

Substituting  equation  (139)  into  (140)  and  the  result  into  (138),  we  have 

r(r)=exp  -^Rdii^i{r)  ,  (141) 


where 

-  ARPzul 

T>  _  ^ 

■th  -  <-) 


(142) 


44 


Hence,  by  plotting  the  MCF  as  a  function  of  the  normalized  range  R  and 
normalized  separation  v  j z,  we  can  generate  universal  sets  of  curves  scal¬ 
able  to  arbitrary  values  of  the  range,  sensor  separation,  frequency,  height, 
friction  velocity,  and  sound  speed. 

The  dependence  of  the  MCFs  on  the  normalized  range,  for  the  Kolmogorov 
model,  various  von  Karman  models,  and  the  Kristensen  et  al  model,  is  plot¬ 
ted  in  figures  15  to  18.  The  first  two  plots  show  along-wind  MCFs  for  sensor 
separations  equal  to  0.1  and  1  times  the  height;  the  second  two  are  the  cor¬ 
responding  crosswmd  plots.  When  r/z  =  0.1  ,  the  various  models  agree 
fairly  well.  The  MCF  for  the  Kolmogorov  model  is  only  slightly  less  than 
the  others.  The  differences  become  more  pronounced  when  r/z  =  1.  At  this 
separation,  the  Kolmogorov  model  predicts  much  lower  coherence  than  the 
others.  Also,  the  isotropic  von  Karman  model  predicts  a  coherence  that  is 
too  high  in  the  along-wind  direction. 

A  similar  set  of  comparisons,  this  time  with  the  Kolmogorov,  Gaussian, 
Kristensen  et  al,  and  Mann  models,  is  shown  in  figures  19  to  22.  When 
rjz  =  0.1,  the  predicted  coherence  from  the  Gaussian  model  is  much  too 
high.  This  is  because  the  Gaussian  model  does  not  have  a  realistic  inertial 
subrange.  The  Gaussian  model  works  well  when  rjz  =  1,  though.  MGFs 
from  the  Kristensen  et  al  and  Mann  models  agree  reasonably  well,  although 
for  crosswind  propagation  when  rjz  =  1,  the  Mann  model  predicts  signif¬ 
icantly  higher  coherence. 


Figure  15.  MCF  as  a  function  of  normalized  range,  for  Kolmogorov,  various  von  Karman, 
and  Kristensen  et  al  models.  Propagation  is  along  wind.  Sensor  separation  is  0.1  times 
height. 
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Figure  16.  MCF  as  a  function  of  normalized  range,  for  Kolmogorov,  various  von  Karman, 
and  Kristensen  et  al  models.  Propagation  is  along  wind.  Sensor  separation  equals  height. 


Figure  17.  MCF  as  a  function  of  normalized  range,  for  Kolmogorov,  various  von  Karman, 
and  Kristensen  et  al  models.  Propagation  is  crosswind.  Sensor  separation  is  0.1  times  height. 


46 


Normalized  range,  R 

Figure  18.  MCF  as  a  function  of  normalized  range,  for  Kolmogorov,  various  von  Karmm 
models,  and  Kristensen  et  al  models.  Propagation  is  crosswind.  Sensor  separation  equals 
height. 


Figure  19.  MCF  as  a  function  of  normalized  range,  for  Kolmogorov,  isotropic  von  Karman, 
Gaussian,  Kristensen  et  al,  and  Mann  models.  Propagation  is  along  wind.  Sensor  separation 
is  0.1  times  height. 
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Figure  20.  MCF  as  a  function  of  normalized  range,  for  Kolmogorov,  isotropic  von  Karman, 
Gaussian,  Kristensen  et  al,  and  Mann  models.  Propagation  is  along  wind.  Sensor  separation 
equals  the  height. 


Figure  21.  MCF  as  a  function  of  normalized  range,  for  Kolmogorov,  isotropic  von  Karman, 
Gaussian,  Kristensen  et  al,  and  Mann  models.  Propagation  is  crosswind.  Sensor  separation 
is  0.1  times  height. 
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Figure  22.  MCF  as  a  function  of  normalized  range,  for  Kolmogorov,  isotropic  von  Karman, 
Gaussian,  Kristensen  et  al,  and  Mann  models.  Propagation  is  crosswind.  Sensor  separation 
equals  height. 


9.2  Performance  of  Acoustic  Direction-Finding  Arrays 

An  important  application  of  the  theory  of  mutual  coherence  is  the  pre¬ 
diction  of  source  direction-finding  accuracy  from  acoustical  arrays.  A  full 
treatment  of  this  problem  was  provided  by  Wilson  (1997a).  The  procedure 
described  in  that  report  involved  calculating  the  Cramer-Rao  lower  bound 
(CRLB)  corresponding  to  a  given  array  configuration,  turbulence  model, 
and  noise  environment.  The  CRLB  is  the  standard  deviation  in  the  wave- 
front  angle-of-arrival  estimates  about  their  actual  value  attained  using  an 
ideal  estimation  method.*  Or  more  simply  put,  when  the  CRLB  is  equal  to 
a  degrees,  we  can  usually  measure  the  direction  of  a  source  to  an  accuracy 
of  a  degrees,  if  equipment  and  processing  methods  are  good. 

Actual  CRLB  calculations,  for  many  of  the  models  discussed  earlier  in  this 
report,  are  shown  in  figures  23  to  30.  The  calculations  were  performed  for  a 
six-element  array,  with  the  sensors  evenly  spaced  around  a  horizontal  cir¬ 
cle  having  a  radius  of  1.2  m.  It  turns  out  that  the  CRLB  for  such  an  array 
configuration  is  independent  of  the  source  direction  (Wilson,  1997a).  The 

‘The  CRLB  is  strictly  the  variance  of  the  estimates,  not  their  standard  deviation  (square 
root  of  the  variance).  Because  the  standard  deviation  has  linear  units  and  is  therefore  more 
natural  to  deal  with,  however,  1  have  adopted  the  implicit  convention  of  meaning  the  square 
root  of  the  CRLB  whenever  1  refer  to  the  CRLB  throughout  this  report. 
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source  and  array  heights  were  both  1  m,  the  friction  velocity  was  tt*  =  0.5 
m/ s  (windy  conditions),  the  SNR  was  20  dB,  and  five  statistically  indepen¬ 
dent  measurements  were  used  for  each  angle-of-arrival  estimate.  The  plots 
show  the  dependence  of  the  CRLB  for  horizontal  bearings  on  source  range 
(propagation  distance)  and  frequency. 

The  general  appearance  of  the  CRLB  plots,  for  all  of  the  models,  is  similar. 
Near  the  source  and  at  low  frequency,  the  CRLB  is  independent  of  range. 
This  independence  is  because  the  main  determinant  of  array  performance 
near  the  source  is  the  environmental  noise.  As  the  frequency  and  range  are 
increased,  however,  the  CRLB  contours  turn  upward  dramatically.  Within 
this  region,  distortion  of  the  wavefronts  by  turbulence  begins  to  be  the  most 
significant  determinant  of  array  performance. 

The  Kolmogorov  model  (fig.  23)  predicts  higher  values  for  the  CRLB  than 
the  other  models,  particularly  at  lower  frequencies  and  longer  ranges.  The 
reason  is  that  the  Kolmogorov  model  predicts  too  much  energy  at  the  large 
scales,  which  are  most  important  in  determining  the  array  performance. 
Therefore  the  Kolmogorov  model  should  be  avoided  in  this  application. 
CRLB  predictions  for  the  isotropic  von  Karman  model  (fig.  24)  agree  rea¬ 
sonably  well  with  calculations  from  the  more  complicated,  anisotropic  mod¬ 
els  (figs.  25  to  30).  Its  predictions  fall  midway  between  the  along-wind 
and  crosswmd  predictions  of  the  Kristensen  et  al  (figs.  26  to  27)  and  Mann 
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Figure  23.  CRLB  calculations  for  Kolmogorov  inertial-subrange  model. 
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Figure  24.  CRLB  calculations  for  isotropic  von  Karman  model. 


Propagation  distance  (m) 

Figure  25.  CRLB  calculations  for  anisotropic  bottom-up  von  Karman  model,  for  propaga¬ 
tion  in  along-wind  direction. 
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Figure  26.  CRLB  calculations  for  Kristensen  et  al  model,  for  propagation  in  along-wind 
direction. 


Figure  27.  CRLB  calculations  for  Krisfensen  et  al  model,  for  horizontal  propagation  in 
crosswind  direction. 
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Figure  28.  CRLB  calculations  for  bottom-up  anisotropic  Gaussian  model,  for  propagation  in 
along- wind  directions. 


Figure  29.  CRLB  calculations  for  Mann  model,  for  propagation  in  along-wind  direction. 
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Figure  30.  CRLB  calculations  for  Mann  model,  for  horizontal  propagation  in  crosswind 
direction. 

(figs.  29  to  30)  models.  Predictions  from  the  bottom-up  von  Karman  model, 
in  the  along-wind  direction  (fig.  25),  are  close  to  the  corresponding  predic¬ 
tions  for  the  Kristensen  et  al  (fig.  26)  and  Mann  (fig.  29)  models.  The  bottom- 
up  Gaussian  model,  in  the  along-wind  direction  (fig.  28),  also  agrees  well 
with  the  corresponding  predictions  from  the  Kristensen  et  al  and  Mann 
models,  although  the  contours  start  to  diverge  for  high  frequencies  and 
long  propagation  distances.  This  may  be  due  to  the  inability  of  the  Gaus¬ 
sian  model  to  realistically  capture  small-scale  turbulence  structure. 

The  Kristensen  et  al  and  Mann  models  are  probably  the  most  realistic,  and 
agree  in  most  regards.  Interestingly,  the  anisotropic  effect  predicted  by  the 
Mann  model  is  more  pronounced  than  the  effect  predicted  by  the  Kris¬ 
tensen  et  al  model.  In  the  along-wind  direction,  the  Mann  model  (fig.  29) 
predicts  higher  GRLB  than  the  Kristensen  et  al  model  (fig.  26);  in  the  cross- 
wind  direction,  it  predicts  lower  GRLB  (figs.  27  and  30). 

On  the  basis  of  the  GRLB  calculations,  we  are  led  to  the  following  conclusions: 

1.  The  Kolmogorov  model  does  not  work  very  well  for  low-frequency 
acoustics  and  should  therefore  be  avoided. 

2.  The  isotropic  von  Karman  model  is  a  convenient  and  reasonably  ac¬ 
curate  model  for  rough  calculations. 
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3.  The  anisotropic  versions  of  the  von  Karman  and  Gaussian  models 
have  little  practical  value.  If  one  is  willing  to  go  through  the  effort 
of  using  an  anisotropic  von  Karman  or  Gaussian  model,  it  would  be 
worthwhile  to  take  the  small  additional  step  to  the  Kristensen  et  al 
model,  or  even  consider  using  the  Mann  model. 
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10.  Concluding  Remarks 


The  assumption  of  isotropic  turbulence  is  a  common  one  and  a  reason¬ 
able  starting  point  for  creating  statistical  turbulence  models.  Previous  re¬ 
searchers  have  used  this  assumption  to  derive  equations  for  the  second- 
order  statistics  needed  in  wave  propagation  calculations.  By  starting  with 
a  Gaussian-type  equation  for  the  correlation  function,  for  example,  or  with 
a  von  Karman-type  equation  for  the  energy  spectrum,  one  can  use  well- 
known  relationships  for  isotropic  homogeneous  turbulence  to  derive  any 
other  second-order  statistics  of  interest. 

The  main  motivation  for  selecting  the  isotropic  Gaussian  and  von  Karman 
models  (as  opposed  to  Kolmogorov's  inertial-subrange  model)  is  that  they 
are  well-behaved  for  the  large-scale,  energy  subrange  turbulence.  This 
large-scale  turbulence  is  important  in  low-frequency  acoustics.  Although 
the  isotropic  Gaussian  and  von  Karman  models  are  more  satisfactory  for 
large-scale  turbulence  than  the  Kolmogorov  model,  they  are  still  unrealistic 
in  one  critical  aspect:  they  neglect  anisotropy  of  the  large-scale  turbulence. 

My  main  intention  in  this  report  was  to  remedy  this  problem,  by  show¬ 
ing  how  anisotropic  versions  of  the  Gaussian  and  von  Karman  models 
can  be  built.  This  task  proved  surprisingly  difficult.  Even  though  a  sim¬ 
ple,  idealized  form  of  anisotropy  was  assumed  (eq  (23)),  the  mathematical 
development  of  the  models  quickly  became  complex.  Difficulties  in  deriv¬ 
ing  closed-form  equations  for  statistical  quantities  of  interest,  such  as  the 
2D  correlation  functions,  made  numerical  integrations  necessary. 

The  inherent  complexity  of  even  "simple"  anisotropic  models  with  obvious 
shortcomings  (such  as  those  developed  based  on  eq  (23))  makes  their  direct 
application  unappealing.  In  practice  it  is  probably  not  much  more  difficult 
to  use  the  first-principle  model  developed  by  Mann  (1994). 

This  report  has  featured  extensive  comparisons  between  different  statistical 
turbulence  models,  with  and  without  anisotropy.  Although  such  compar¬ 
isons  are  a  helpful  first  step  toward  understanding  the  effects  of  anisotropic 
turbulence  on  acoustic  propagation,  in  the  end  experiments  are  needed  to 
test  how  well  the  various  anisotropic  models  actually  work  for  predicting 
acoustic  propagation  characteristics.  Such  experiments  are  left  for  future 
research. 
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Appendix  A.  Some  Useful  Integrals 


Several  solutions  to  integrals  that  were  used  in  the  body  of  this  report  are 
provided  in  this  appendix,  along  with  citations  of  their  sources. 

1.  From  Beyer  (1981),  definite  integral  equation  (615), 


rcxi  y,n  1^^  _  r  (n/2)  r  (v  —  n/2) 
'o  (l  +  x2)"  =  2T{u) 


]^B{n/2,u-n/2).  (A-1) 


roo  y2m+l 

Jo  {c  +  ay’^Y 


By  changing  the  variable  of  integration  from  x  to  y  =  cx/a  in  the 
above  equation  and  setting  n  =  2m  +  2,  one  arrives  at  the  useful 
equation 


r(m  +  l)r(;y-m-l)  ^  1  R  +  1  _  rr.  -  1 1 

2r  (v)  ’ 


(A-2) 


2.  From  Gradshteyn  and  Ryzhik  (1994),  equation  (6.565.4), 


roo 


3.  From  Gradshteyn  and  Ryzhik,  equation  (3.194.1), 

rs  3  ^iy-\-£—2  /  s\ 


where  2F1  is  the  hypergeometric  function.  Actually,  the  integral  can 
be  rewritten  in  terms  of  an  incomplete  beta  function  rather  than 
the  hypergeometric  function.  This  is  convenient,  since  many  soft¬ 
ware  packages  such  as  Matlab  have  routines  for  calculating  incom¬ 
plete  beta  functions,  but  not  hypergeometric  functions.  It  should  be 
pointed  out  that  these  routines  usually  require  that  the  argument  to 
the  incomplete  beta  function  is  in  the  range  [0, 1];  hence  we  must  try 
to  recast  our  result  in  a  form  satisfying  this  requirement.  Let  us  begin 
by  substituting  u  =  into  the  right-hand  side  of  equation  (A-4). 
We  find 


h  (i/<2 +  !)'■'“  k 


■s/p  yU+e-'i 

,  {u+lf 
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Next  we  adopt  a  trick  from  Kristensen  et  al  (1989)  by  setting  tj  = 
u/  {1  +  u),  with  result 


Jo  (f/p  +  if 


s/£2/(i+s/£2 


3 


From  the  definition  of  the  incomplete  beta  function  (eq  (84)),  we  have 

PS  fiy+e—S 

Jq  [ijp.  ~  ^ (z^  +  e  —  2, 2  —  e) .  (A- 7) 

Furthermore,  setting  s  =  k~‘^,  we  can  write  this  as 

rk-^  +u+e-‘i 

Jq  [ijp.  ~  {u  +  e  —  2,2  —  e) .  (A-8) 

Note  that  for  all  real  values  of  k,  the  argument  to  the  beta  function, 
1/  (l  +  P£^),  is  indeed  in  the  range  [0, 1]. 

4.  From  Gradshteyn  and  Ryzhik,  equation  (3.381.3), 


p  "T  (v,  fiu)  , 


where  F  (x,  y)  is  the  incomplete  gamma  function.  When  we  substi¬ 
tute,  t  =  1  /x,  and  s  =  1/m,  the  integral  becomes 


5.  In  order  to  calculate  the  2D  correlation  function  for  the  bottom-up 
von  Karman  model,  one  must  solve  integrals  of  the  form 

I  =  (kH^)  (a,  h)  Jo  {kr)  k  dk.  (A-10) 

It  happens  that  these  integrals  can  be  determined  in  terms  of  Mei- 
jer 's  G-functions.  We  begin  by  writing  the  incomplete  beta  function  as 
a  hypergeometric  function  using  equation  (6.6.8)  from  Abramowitz 
and  Stegun  (1965): 

-®i/(i+fc2£2)  (a,  6)  =  a  ^2^2  )  F  ^a,  1  —  6;  a  -I-  1;  ^  ^  ' 
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Next,  we  transform  the  hypergeometric  function  using  equation 
(15.3.4)  from  Abramowitz  and  Stegun: 


-®i/(i+fc2£2)  (a,  6)  —  a  ^  F  ^a,  a  +  6;  a  +  1;  — . 


The  hypergeometric  function  can  now  be  recast  as  a  Meijer's  G- 
function,  using  equations  (9.34.7)  and  (9.31.2)  from  Gradshteyn  and 
Ryzhik  (1994).  The  result  is 


-Si/(1+A:2£2)  (a,  b)  — 


(A:2£2) 


— a— 1 


T  (a  +  6) 


-Gi  I  A:2"2 


2,  a  +  2 

dH-l, 


(A-ll) 


We  can  now  substitute  equation  (A-ll)  into  (A-10)  and  solve  the  in¬ 
tegral  using  equation  (7.822.1)  from  Gradshteyn  and  Ryzhik  (1994). 
The  result  is 


/  = 


2£2r(a  +  6)  \2£) 
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(A-12) 
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